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Summary 

Mixed convection above a horizontal disk rotating in a semi-infinite fluid is examined 
when the disk is heated so that its temperature varies quadratically with distance 
away from its centre. Steady similarity solutions are presented for a range of values 
of the two dimensionless parameters, a Prandtl number and a Grashof number. The 
results are corroborated by asymptotic analyses undertaken in the limits of small and 
large Prandtl number. For finite Prandtl number, the existence of multiple solutions 
at fixed Grashof number is revealed. The similarity structure for steady solutions 
fails beyond a critical Grashof number and this is interpreted in terms of a finite-time 
singularity of the unsteady form of the governing equations. 

1. Introduction 

Free convection refers to the motion created by thermal effects in a fluid which is otherwise 
at rest. Forced convection refers to the situation in which an external flow field, such 
as a uniform stream, directly controls the movement and distribution of heat within a 
fluid. Of particular interest in the present work is the case of mixed convection in which 
free-convective effects interact directly with a forced convective motion. 

In a typical free convection problem, buoyancy forces drive motion directly. This occurs, 
for example, when a heated plate is held inside a fluid in a vertical position. However, 
fluid motion may result from thermal buoyancy effects even when no component of the 
buoyancy force acts in the direction of the subsequent motion. Stewartson (1) and Gill et 
al. (2) demonstrated that free horizontal convection may occur in a boundary layer on a 
semi-infinite horizontal plate at a constant temperature. In this case there is no component 
of the buoyancy force in the direction of motion; rather the horizontal flow occurs indirectly 
as a result of a pressure gradient induced along the plate. Amin & Riley (3) demonstrated 
that horizontal free convection takes place above an infinite plate held in a quiescent fluid 
when the plate temperature varies along its length. This flow is also of boundary-layer type 
provided that the Grashof number, which measures the relative strength of the buoyancy- 
induced convection to viscous effects, is large. 

In a study of mixed convection at a vertical plate held at a constant temperature, Mcrkin 
(4) showed that when the buoyancy force opposes the prevailing direction of motion, the 
boundary layer eventually separates and the solution fails at a singularity (akin to the 
classical Goldstein singularity in an isothermal boundary layer) where the skin friction 
vanishes. Similar singular behaviour and flow separation is seen for uniform flow past a flat 
plate with constant heat flux (5, 6), and a discussion of the character of the singularity 
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at the separation point was provided by Hunt & Wilks (7). Jones (8) showed that for a 
heated semi-infinite plate inclined at a small negative angle to the horizontal, flow separation 
occurs without meeting a singularity, and it is possible to continue the integration of the 
boundary-layer equations downstream. The absence of a singularity may be attributed 
to the ability of the pressure-gradient to adjust interactively to accommodate the flow 
separation in a manner reminiscent of the interactive behaviour of the pressure-gradient 
required in traditional triple-deck theory (e.g. 9). 

A number of mixed convection scenarios have been discussed in which a flow due to 
thermally-induced pressure gradient achieves a balance with a background flow field without 
provoking flow separation. Typically these have been described using similarity structures of 
the stagnation-point type which share many features with traditional boundary-layer theory 
but which do not formally require approximation on the assumption of a small dimcnsionlcss 
parameter. Amin & Riley (10) found that a steady solution may be constructed for mixed 
convection at a stagnation-point on a plane wall when the wall temperature gradient is 
assumed to vary quadratically with distance along the wall. Surprisingly a steady solution 
is possible not only when thermally-driven flow component is in the same direction as 
the background flow but also when it opposes the background flow. Arrieta-Sanagustm 
(11) considered the axisymmetric version of this problem and demonstrated that a steady 
solution is possible with the thermal flow opposing the stagnation-point flow provided that 
the Grashof number lies below a threshold value. Above this value the assumed similarity 
structure fails and an eruption of fluid from the vorticity layer adjacent to the wall is 
anticipated. Further examples are given by Revnic et al. (12), who examined mixed 
convection near a stagnation point on a vertical cylinder, and Merkin & Pop (13) who 
studied mixed convection at a vertical surface exposed to a far-field shear flow. In both cases 
steady solutions are available when the thermal flow reinforces the shear flow. Solutions 
also exist for an opposing thermal flow if a dimensionless mixed convection parameter is 
below a critical value. Finally we mention the very recent work of Riley (14), who examines 
the unsteady flow above a horizontal wall induced by a standing temperature wave. 

In this paper we discuss mixed convection above a heated rotating disk. In the absence 
of a temperature field, the flow is described by a set of ordinary differential equations 
according to the classical three-dimensional stagnation-point structure of von Karman (e.g. 
15). Solutions of these equations have been studied exhaustively. Cochran (16) computed 
the most well-known solution, and that which is most commonly described in textbooks 
(e.g. 15, 17). More recently, a number of authors have pointed out the existence of multiple 
solutions. Zandbcrgcn and Dijkstra (18) computed a second solution as part of a wider 
study which accounted for the effect of suction at the plate. Later, Lentini and Keller (19) 
showed that in fact there are an infinite number of solutions. An excellent survey of the 
results is given by Zandbergen and Dijkstra (20). In the presence of a temperature field, the 
flow can be analysed using the same von Karman similarity structure with the additional 
assumption that the wall temperature varies quadratically with distance away from the 
centre of the disk. The mixed convection is characterised by two dimensionless parameters, 
a Prandtl number and a Grashof number. The rotation of the disk draws fluid downwards 
from infinity and forces it outwards away from the disk centre. For positive Grashof number, 
the thermally-induced pressure gradient opposes this and attempts to drive fluid towards 
the axis of rotation. For negative Grashof number, the thermally-induced pressure gradient 
complements the centrifugal flow. This case was studied by Sreenivasan (21). 
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Section 2 is devoted to steady solutions for a fixed Prandtl number, which are expected 
to exist when the Grashof number lies below a critical value. Moreover the non-uniqueness 
mentioned above for isothermal flow strongly suggests the presence of multiple solution 
branches for the heated flow problem, and this is indeed found to be the case. The steady 
solutions are investigated further in the limit of small and large Prandtl number in section 
3. Above the critical Grashof number the solution to the unsteady, self-similar form of the 
governing equations encounters a finite time singularity, and this is explored numerically in 
section 4. Our findings are summarised in section 5. 

2. Steady flow 

We consider the steady flow generated by a combination of a horizontal, rotating flat plate 
and a thermally-induced radial pressure gradient. The flow is assumed to be incompressible 
and satisfy the Navier-Stokes equations under the Boussinesq approximation. The rotating 
plate boundary is located at z* = and the fluid occupies the region z* > 0. The plate is 
assumed to rotate about the z* axis with angular velocity ilk, where k is the unit vector in 
the positive z* direction. The flow is assumed to be axisymmetric with velocity components 
u*,w* in the r* , z* directions respectively and v* is the azimuthal velocity component. The 
pressure in the fluid is represented by p* and the temperature field in the fluid is represented 
by 9* . The plate temperature is assumed to vary quadratically with distance from the origin 
so that 

e* = e* -( K ^)(o*o-oio)r 2 (2.1) 

at z* = 0, where b is a positive dimensionless parameter, 9*^ is the ambient temperature 
far from the plate and 9q is the plate temperature at the origin. 
We non-dimensionalize by writing 

(r*,z*) = (u/n) 1 ^(r,z), (u*,v*,w*) = (M) 1/2 (u,v,w), (2.2) 
p* = (pM)p, e* = 0* oo + [61-6^)6, (2.3) 

where v is the kinematic viscosity. According to the Boussinesq approximation the fluid 
density is given by 

P = Poo(i-P(9* -e* oo )e), (2.4) 

where p^ is the density attained at the ambient temperature 9*^ and j3 is the coefficient of 
thermal expansion far from the plate. The temperature field inside the fluid is governed by 
the dimensionless energy equation 

cr(u9 r + w9 z ) =9 rr + 9 r /r + 9 zz , (2.5) 

where a — v / n is the Prandtl number and k is the thermal diffusivity of the fluid. 
We seek a similarity solution of the form 



u = ru (z), v = rv (z), w = w (z), p = p (z)+r 2 p 1 (z), 9 = 9 (z) + br 2 9 1 (z). (2.6) 
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Substituting (2.6) into the Navier-Stokes equations we obtain 



2u + ^° = 0, (2.7) 
dz 

w o +w o^=-2pi + ^, (2.8) 

d^p _ djvo 
dz dz 2 



2u v + w — = (2.9) 



2uo e 1 + Wo ^i = ^. (2 . 10) 

dz a dz z 

The boundary conditions are 

u = 0, v = 1, w = 0, 0i = -1 at z = 0, (2.11) 

and 

Mo, ^o, 0i — > as z — > oo. (2.12) 
The pressure field pi is given by 

/•OO 

pi = -Gr / 0i dz, (2.13) 

J z 

where the dimensionless Grashof number is defined to be 

Gr ~ ^i/2 3/2 • ( 2J4 ) 

The remaining unknown functions in (2.6), namely Po(z) and 0o(z), satisfy a pair of ordinary 
differential equations, with accompanying boundary conditions, which can be written down 
in a straightforward manner. These arc not the main focus of the present work, and so we 
do not give these here; our primary concern is with the system (2.7)-(2.13). 

The problem is solved numerically by first truncating the infinite domain of integration 
at a finite level, z^, and introducing a non-uniform grid covering the range [0, Zoo]- In 
pratice it is useful to cluster grid points at z = in order to more accurately resolve the 
solution in the region close to the wall where variation is typically rapid. The equations are 
discretized using centered differences for the derivative terms at each interior grid point. 
The boundary conditions (2.11) are enforced at the first grid point and the conditions in 
(2.12) together with pi = are enforced at the last grid point. The non-linear terms in 
the governing equations are dealt with using a quasi-linearization procedure as part of an 
iterative process. For example the term Uq in (2.8) is replaced by 2uqUq — u 2 , where uq is 
the value from the previous iteration. In this way at each iteration we solve a tri-diagonal 
system of linear algebraic equations for the unknowns at each grid point. This is done 
efficiently using the Thomas algorithm (e.g. 22). The solution is deemed to have converged 
once the condition |uo — Uo\ < S is met at every grid point for a prescribed choice of S. A 
similar convergence condition is required for all of the other flow variables. In practice we 
took S = 1CT 8 . As a check, the converged numerical solution is substituted back into the 
governing equations. 
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Fig. 1: The steady solution curves for a = 1. The main branch is shown with a solid line, 
and one of the secondary branches is shown with a broken line. Velocity and temperature 
profiles at the points labelled 1, 2, 3 are shown in figures 2, 3 and 4 respectively. 



When Gr = 0, so that thermal effects are absent, the problem reduces to the classical 
von Karman rotating disk flow for which fluid is swept radially outwards along the plate 
and drawn down from infinity along the axis of rotation. For Gr > a thermally-induced 
pressure gradient tends to force fluid radially inwards and up away from the plate. The 
presently sought similarity solution assumes that the magnitude of the vorticity decays 
with vertical distance upwards from the boundary. Its existence will be determined by 
a competition between the centrifugally-driven downwash, which tends to draw vorticity 
towards the boundary, and the thermally-driven upwash which tends to carry it away from 
the boundary. If Gr is sufficiently large we expect the latter to overcome the former and, 
consequently, for the hypothesised similarity structure to fail. 

As was mentioned in the Introduction, the presence of multiple solution branches is 
expected. In Fig. 1 we show two of the solution branches for the case a = 1. The flow is 
characterised using the negative of the vertical velocity at infinity, Woo. On both branches 
the solution at Gr = was computed using the numerical method described above. The 
branches were then extended to positive Grashof numbers using the parameter continuation 
software AUTO-07p (23) with the finite difference solution at Gr = used as a starting 
point. The branch shown with a solid line starts at the point (Gr, Woo) = (0,0.884), which 
coincides with the appropriate calculation of Zandbergen and Dijkstra (18) (see their table 
4.1). The branch shown with a broken line starts at (Gr, w^) = (0,0.447723), which also 
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agrees with table 4.1 of Zandbergen and Dijkstra (18). As the branch extends away from 
Gr = and starts to oscillate down toward the Gr axis, we found that it was necessary to 
gradually increase as the thickness of the vorticity layer at the disk grows. Over most 
of the branch we took = 500, but towards the lowermost part of the branch shown in 
the figure it was necessary to take the very large value Zoo — 2 x 10 4 . For practical reasons, 
the computations were terminated at the end-point shown. It seems reasonable to suppose 
that the branch will continue to approach the Gr axis while the thickness of the vorticity 
layer continues to grow. 

Having delineated the main branch as described, a few individual points along the branch 
were confirmed independently using an alternative numerical method as follows. The value 
of Woo was fixed and Gr was treated as an unknown to come as part of the solution. 
Starting from the large z asymptotic form of the solution provided by (18), (2.7)-(2.10) were 
integrated inwards from infinity, and a shooting method was used to satisfy the conditions 
at the disk (2.11). In each case the resulting pair (Gr, — Woo) was confirmed to lie on the 
main branch in Fig. 1. 

On the main solution branch shown in Fig. 1, there is one solution when < Gr < 0.0156, 
there are three solutions when Gr > 0.0156, and there are no solutions when Gr > Gr c , 
where the critical Grashof number Gr c = 0.0466. Velocity and temperature profiles at 
three locations on the main branch are shown in Figs 2, 3 and 4. Note that in the presently 
considered case of a — 1, it is readily seen from the governing equations (2.9) and (2.10) 
and the boundary conditions (2.11) and (2.12) that vq = —9\ and so the profile for the 
swirl component of velocity may be easily inferred from the temperature profiles shown in 
the figures. In each case the vertical velocity is single-signed and the fluid descends from 
infinity toward the disk. Notably, the radial velocity switches sign in the profiles shown in 
Figs 3 and 4 as a turning point develops in the axial wq profile. In Figs 2, 3, and 4, 6\, and 
hence also the swirl component v a , has the same sign everywhere. 

Other values of the Prandtl number, er, were also investigated, and broadly the same 
behaviour was found as that discussed above. For example, the main solution branch was 
computed for the case a = 0.5 and it is found to be have qualitatively the same shape as 
that shown in Fig. 1. Fig. 5 shows how the critical Grashof number beyond which there 
are no steady solutions depends on the Prandtl number, a. Notably, as a increases, Gr c 
appears to grow without bound, and as a decreases, Gr c appears to approach zero. These 
trends are investigated further in the next section where we examine the limits of large and 
small a in turn. 

3. Steady flow for large and small Prandtl number 

Here we use asymptotic analysis to discuss the flow when a ^ 1 and a <C 1. Of particular 
interest is the behaviour of the the critical Grashof number, Gr c , in these limits. 

3.1 The case a > 1 

For large Prandtl number the diffusivity of heat from the boundary is much less than 
the diffusivity of vorticity created at the boundary. A consequence of this is that the 
thermal boundary layer is much thinner than the velocity boundary layer and so within the 
inner, thermal boundary layer the radial velocity is a linear function of z, i.e. uo cc z. A 
consequence is that if convection and diffusion of heat are to balance therein, (2.10) shows 
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Fig. 2: The velocity profiles wo (solid line), uq (dot-dashed line), and the temperature 
profile Q\ (dashed line) for a — 1 and Gr = 0.0269, corresponding to label 1 in figure 1 with 
Wrx, = —0.803. Note that for u = 1, as here, v = —6\. 



that the thickness of this inner layer is 0(a 3 ), a classic result. Further uq = 0(a 3), 
and the continuity equation (2.7) requires wq = 0(er~ 5 ), with vo = 1 at leading order. A 
balance between the pressure and viscous forces in the radial momentum equation (2.8) 
then requires p\ = 0{a^). Finally we note from equation (2.13) that Grashof numbers as 
large as 0(cr3) are available for steady flow in this situation. For this inner region then we 
introduce the scaled variables 

u a = a~^U , w = a~%W , pi=a*Pi, Gr = cr^A, z = a~^Z, (3.1) 

which yield the following equations for the inner layer 

2(/„ + ^ = 0, = -2P, + gf, (3.2) 

with boundary conditions 

U = W = 0, 6»i = -1 at Z = 0, Pi, <9i as Z -S- 00, (3.4) 



together with a suitable matching condition for Uq as Z — > 00. 
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Fig. 3: The velocity profiles wo (solid line), u (dot-dashed line), and the temperature 
profile 6\ (dashed line) for a = 1 and Gr = 0.0219, corresponding to label 2 in figure 1 with 
Woo = —0.284. Note that for u = 1, as here, vq = —0\. 



Outside this thin thermal boundary layer pi = 6\ = and the velocity components 
are each of order unity. Therefore, the governing equations for the outer layer, whose 
characteristic thickness is of order unity, are from (2.7) to (2.9) 



2u + 



together with 



dw 
~dz~ 



= 0, 



2 2 , dw 

u -v + w a~^ 



d 2 Ug 



2u a v + w —— = 
dz 



d 2 v 
dz 2 ' 



(3.5) 



u = v 



1 



wo 



at z = 0, 



u ,v 



as 



(3.6) 



Therefore the flow in this outer layer corresponds to the classical von Karman rotating disk 
flow. In order to match the inner and outer solutions note from (3.2) that as Z — > oo, 
Uq oc Z. Accordingly, taking account of the scalings in (3.1), the matching requirement is 



dU 
dZ 



duo 

~dz~ 



(3.7) 



z=0 



The system of equations (3.2) to (3.4) has been integrated numerically using the method 
outlined in section 2 for varying values of A. In particular we find that no solutions exist 
for A = A c > 0.436. Therefore the asymptotic behaviour of the critical Grashof number for 
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Fig. 4: The velocity profiles Wo (solid line), uq (dot-dashed line), and the temperature 
profile 6\ (dashed line) for a — 1 and Gr = 0.0176, corresponding to label 3 in figure 1 with 
Woo = —0.122. Note that for u = 1, as here, v = —6i- 



large values of the Prandtl number has 

Gr r - 0.436er3 as a 



oo. 



(3.8) 



This result is compared with the results obtained from the numerical solutions of (2.7) to 
(2.13) in Fig. 5, showing good agreement for large values of a. 



3.2 The case a < 1 

For very small values of the Prandtl number, the thermal diffusivity is much greater than the 
diffusivity of vorticity, with the consequence that the thermal boundary layer now exceeds, 
by far, the velocity boundary layer in thickness. Assuming, as is reasonable, that the range 
of available values of Gr is small (see Fig. 5) in the inner region, where z — 0(1), we have 
from (2.7)-(2.10), at leading order 



dwn 

2u + —^= 0, 
az 



. 2 ~2 , ~ d "0 

u Q -v +™ — 



d 2 ii 



2u a v + w — = -p- , (3.9) 



dz 2 ' 
together with 

= v — 1 = wo = at z = 0, uo, So — > (3.10) 

as z — > oo. Notice that we have used a tilde superscript to indicate variables in the inner 
region. System (3.9), (3.10) is uncoupled from the temperature field and is identical to 
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Fig. 5: The dependence of the critical Grashof number, Gr c , on the Prandtl number a, 
shown with a solid line. The broken lines correspond to the asymptotic predictions, Gr c = 
0.436cr 2 / 3 in (3.8) and Gr c = 0.0123er 3 in (3.24), for large and small a respectively. 



that which governs classical von Karman flow over a rotating disk. The solution has the 
property that w n — » — k as z — » oo, where k = 0.884 (e.g. (20), table 4.1), and u , «o — > 
exponentially fast. 

In the outer thermal boundary layer, we expect a balance between the outward diffusion 
and the inward convection of heat so that in (2.10) 

ctw o-t- = \ak— . (3.11) 



dz 2 dz \ dz 

This suggests that the aforementioned balance occurs in a region of scale 0(a~ 1 ) which 
prompts the definition of an outer variable ( = az. Written in terms of the new outer 
variable, (2.8) and (2.13) become respectively, 

2 2 , du o „ . 2 d2 "o dpi Gr 
u -v -ak— = -2 Pl +a = ~ (3-12) 



The first of these suggests viscous diffusion is unimportant and, balancing the three 
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momentum terms on the left hand side with the pressure term on the right hand side, 
indicates that u = 0(a), vo = 0(a) and p\ = 0(a 2 ) in this outer region. Substituting the 
last of these scalings into the second equation in (3.12), and assuming that 9\ = 0(1), we 
deduce that Gr = 0(a 3 ). Accordingly we now write Gr = Act 3 where A = O(l) and seek to 
determine the critical value of A beyond which there is no steady solution. 

With the above preliminaries in place we now develop the solution in greater detail. In 
the inner region we expand the velocity, pressure and temperature fields by writing 

u a = u + aui + a 2 u 2 H , v = % + avi + a 2 v 2 H , 

w = w + awi + a 2 w 2 H , pi — a 2 p H , #i = #i + a9 2 H . (3.13) 

Introducing these expansions into (2.7)-(2.10) we obtain (3.9) at leading order for which, 
as was mentioned above, the classical von Karman solution is available. The leading order 
temperature field satisfies 



d 2 <9i 
~dz^ 



= 0, (3.14) 



with 9\ = — 1 on z = 0. Integrating and applying the wall boundary condition we find 
9i = —1 + Az for constant A. We find that a self-consistent solution can be constructed by 
taking A = so that 9\ = —1. At first order, namely O(a), we obtain from (2.7) and (2.8) 
respectively, 

dtth „_ _ _ _ dui _ duo d 2 ui 

2ui + — — = 0, 2u mi - 2w «i + w — h = -j -s-, (3.15) 

dz dz dz dz A 

and from (2.9) we have 

n ~ ~ n ~ ~ ~ dwi . dv d 2 vi , n , „. 

2u «i + 2miw + wo— hwi— = ——. (3.16) 

dz dz dz z 

At the wall we have u\ = v\ = w\ = 0. Taking the limit z — > oo, and considering the 
dominant balances of terms in (3.15), (3.16) we find 

u\ <~ ai + e.d.t. wi ~ — 2ctiz + a 2 + c.d.t., v\ ~ + c.d.t. (3.17) 

as z — > oo, where e.d.t. means 'exponentially decaying terms' and a\, a 2 , and are 
constants. 

In the outer region, we expand by writing 

u = a &(() + ■■■ , Vo = aV(0 + --- , w = W(C) + ---, (3.18) 
Pl =a 2 P(C) + --- , 1 = e(C) + ---, (3.19) 
where C = crz. Substituting into (2.7)-(2.10) and (2.13) we obtain at leading order, 

2U+ d ^=0, U 2 - V 2 + = -2P, (3.20) 

+ f = A6, « + *|-« (3.21) 



12 J. Arrieta-Sanagusti'n, M. G. Blyth 




Fig. 6: ij<1: The solution branches for the first order problem in the inner region (lower 
panel, showing r = dwi/dz| z= o) and the leading order problem in the outer region (upper 
panel, showing U(0)). 



together with 

W(0) = -k, 0(0) = -1, and U,V,P,O^0 as C -> oo. (3.22) 

The first two conditions in (3.22) are required to match to the inner region. Moreover, 
matching is completed by taking 

ai - U(0), a 3 = V(0). (3.23) 

Integrating the first equation in (3.21), with the use of the first equation in (3.20), we 
find that V = — (a 3 /k)W, where the constant of integration has been chosen to fulfil the 
matching condition in (3.23). Substituting this result, the system (3.20)-(3.22) is solved 
numerically using a shooting method based on fourth-order Runge-Kutta integration. The 
numerical solution shows that az — so that V = 0. The numerical solution is illustrated in 
the upper panel of Fig. 6, where U(0) is plotted against A. Evidently two possible solutions 
exist over the range of A values shown, and there is no solution beyond the critical value 
A c = 0.0123. The lower panel of the same figure shows the corresponding solution branch, 
showing dui/dz\ z = versus A, for the first order inner system (3.15)-(3.17). It should be 
noted that while the solution for the swirl component V is identically zero in the upper 
region, the inner solution V\ varies smoothly from zero at z = 0, through non-zero values, 
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toward zero as z — > oo. According to the preceding remarks, we now have the asymptotic 
estimate for the critical Grashof number, 

Gr c - 0.0123cr 3 as a -> 0. (3.24) 

This is compared with the critical value of the Grashof number obtained from the full 
equations in Fig. 6, showing excellent agreement. 

In concluding this section it is worth reflecting on the rather unexpected flow structure 
which has been uncovered. First, we highlight the somewhat surprising form of the inner 
layer expansions (3.13). As was noted above, the leading order velocity field, namely 
(uo, v , wo), coincides with that for classical von Karman flow on a rotating disk. Deviation 
from this basic profile occurs as a result of a thermally-induced pressure gradient which 
arises from the spatially-dependent wall temperature profile. It is notable that although this 
thermally-induced pressure gradient appears at 0{a 2 ) in the inner region expansions (3.13), 
it is necessary to include velocity perturbations of size 0(a). These velocity perturbations 
are driven indirectly by thermal effects via the radial slip velocity which is imposed upon 
the inner layer by the solution in the outer region, where the effect of a thermally-induced 
pressure-gradient is felt at leading order. Finally we note that the present flow structure 
is similar in vein to that uncovered by Riley (14) in his study of the flow response of a 
semi-infinite fluid above a horizontal wall with a standing wave temperature profile. 



4. Time-dependent calculations 

Following the remarks of the previous section, it appears that no steady solution exists 
when the Grashof number lies above a critical value, Gr c , which depends on the size of 
the Prandtl number. To explore the flow behaviour above the critical Grashof number, we 
consider the unsteady form of the governing equations (2.7)-(2.13) with the time-dependent 
terms included, 

2u + ^ = 0, (4.1) 

^+^-,o 2 +-o^ = -2, 1 + ^#, (4.2) 
9vq dv a d 2 v 

— +2u oVo + Wo — = —, (4.3) 

90i n „ &e 1 i d 2 e 1 

-+2^+^- = -^, (4.4) 

where t = Qt* is a dimensionless time, together with (2.13). These are to be integrated 
forwards in time subject to the boundary conditions (2.11) and (2.12). A suitable initial 
condition is given by 

1 dF 

U °"~2cb' w o = G, w = F, 1= p 1= O (4.5) 

at t = 0, where the functions F(z) and G(z) are Cochran's (16) solution to (2.7)-(2.9) at 
Gr = 0. We note that the same initial condition was used by Riley (24, 25) and Andrews 
and Riley (26) in their studies of heat transfer on a rotating disk. 
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Fig. 7: Velocity and temperature profiles u , 
for cr = 1 and Gr = 0.05. 




vq, wo, and 6\ at various dimensionless times 



The integration forwards in time may be achieved numerically employing a Crank- 
Nicolson scheme based on the same discretisation of the equations as in Bodonyi and 
Stewartson (27). The nonlinear terms are dealt with at each time step using an iterative 
method similar to that described by Tanehill et al. (28). A uniform time-step of size 
0.01 was used for the Crank-Nicolson scheme, whereas the same spatial discretization as 
described in section 2 was employed. 

As in the previous section we focus attention on the case a = 1. The velocity and 
temperature profiles uq, vo, wq and 9\ are shown in Fig. 7 at different times t. Similar 
behaviour is encountered at other Grashof numbers above the critical value. In the early 
part of the flow evolution, the majority of the vorticity is confined to a layer adjacent to the 
disk. As time increases, the thickness of this layer increases as vorticity spreads further and 
further from the disk; this is evident from the behaviour of the v curves in Fig. 7, which 
approach zero further and further from the disk as time goes on. Eventually, the numerical 
evidence is that a finite-time singularity is encountered and the solution blows up (in this 
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case at t w 107). This singular behaviour lends strong support to the assertion that no 
steady solution exists above the critical value of the Grashof number discussed in section 2. 

5. Summary 

We have studied mixed convection over a rotating disk with the flow driven in part by 
the rotation of the disk and in part by a thermally-induced pressure gradient which arises 
in response to a prescribed temperature profile in the disk. Working on the basis of the 
Boussinesq approximation, we derived a set of partial differential equations describing the 
flow and the temperature profile throughout the fluid according to an assumed similarity 
structure. Steady solutions of this system are possible when the dimensionless Grashof 
number lies below a critical value which depends on the size of the Prandtl number. 
Below the critical Grashof number, an infinite number of solution branches exist, each 
one stemming from one of the infinity of branches which exist at zero Grashof number for 
isothermal flow. The main solution branch, which stems from the well-known Cochran 
solution of the classical von Karman problem at zero Grashof number, undergoes a 
number of successive turns so that multiple steady solutions exist along this branch alone. 
Asymptotic descriptions of the flow, and in particular estimates for the critical Grashof 
number, were obtained in the limit of large and small Prandtl numbers and shown to 
be in agreement with the numerical computations. Beyond the critical Grashof number 
there are no steady solutions; rather, our unsteady calculations suggest that the solution 
of the unsteady similarity equations, subject to a suitable initial condition, terminates in a 
finite-time singularity. 
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